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The problem of infrared divergence of tlie effective electromagnetic field produced by elementary 
charges is revisited using the model of an electron freely evolving in a photon bath. It is shown 
that for any finite travel time, the effective field of the electron is infrared-finite, and that in each 
order of perturbation theory the radiative contributions grow without bound in the large-time 
limit. Using the Schwinger-Keldysh formalism, factorization of divergent contributions in multi-loop 
diagrams is proved, and summation of the resulting infinite series is performed. It is demonstrated 
that the effective electromagnetic field of the electron vanishes in the large-time limit, and that 
this vanishing respects the total charge conservation and the Gauss law. It is concluded that the 
physical meaning of infrared singularity in the effective field is the existence of a peculiar irreversible 
spreading of electric charges, caused by their interaction with the photon field. This spreading 
exists in vacuum as well as at finite temperature, and shows itself in a damping of the off-diagonal 
elements of the momentum-space density matrix of electron, but does not affect its momentum 
probability distribution. It precludes preparation of spatially localized particle states at finite times 
by operating with free particle states in the remote past. Relationship of the obtained results to the 
Bloch-Nordsieck theorem is established, and discussed from the standpoint of measurability of the 
electromagnetic field. The effect of irreversible spreading on the electron diffraction in the classic 
two-slit experiment is determined, and is shown to be detectable in principle by modern devices 
already at room temperature. 

PACS numbers: 12.20.-m, ll.lO.Wx, ll.15.Bt, 03.65.Yz 

I. INTRODUCTION 

The low-energy behavior of quantum systems is an issue which is important in both fundamental and applied aspects 
of quantum field theory. Investigation of the low-energy properties of particle interactions constitutes an essential 
part in establishing the correspondence between classical and quantum theories, and finds numerous applications in 
all areas of quantum physics, from noise theory to quantum cosmology. As is well known, these properties essentially 
depend on the presence of massless particle states. In non-relativistic quantum mechanics, the existence of massless 
photons and gravitons shows itself in the form of long-range Coulomb and Newton forces between massive particles. 
Account of the radiation effects brings in, among other outcomes of the union of quantum theory with the relativistic 
principle, a factor that has no counterpart in the physics of systems involving only massive particles. Namely, 
the masslessness of photons and gravitons, together with the non-conservation of the particle number, imply that 
production of these quanta in any scattering process is beyond the experimental control. More precisely, finite 
sensitivity of any experimental setup does not allow one to distinguish scattering processes which involve different 
numbers of massless quanta with sufficiently small energy. 

The role of this indistinguishability in the scattering theory is also known well. In the standard formulation using 
the S-matrix, the scattering process is considered formally on an infinite time interval, implying that the 4-momenta of 
virtual photons describing radiative corrections may take on arbitrarily small values. Integration over such momenta 
gives formally infinite results at every order of the perturbation theory. On the other hand, the uncontrollable 
production of photons with arbitrarily small energy means that the observed scattering cross-section is actually an 
infinite sum of terms, each of which represents the given scattering process with a definite number of extra real soft 
photons, i.e., photons with energies below the sensitivity threshold.^ Integration of the cross-section over the real soft 
photon momenta brings in another divergence which exactly cancels the divergence due to virtual photons, resolving 
thereby the infrared catastrophe of quantum electrodynamics |1| . For brevity, this result will be referred to below as 
the Bloch-Nordsieck theorem. Similar cancelations take place in quantum gravity , and in a more intricate way, in 
Yang-Mills theories ^ . 

These results, though resolve the infrared catastrophe in a physically adequate way, leave open the question about 
possible physical manifestations of the infrared singularities, other than mere explicit dependence of the cross-sections 



And with a total energy going into unobserved photons loss than the uncertainty in the energy of "hard" particles, i.e., the particles 
being scattered. 
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on the sensitivity threshold. That the Bloch-Nordsieck theorem does not exhaust the infrared problem can be seen 
from the standpoint of the measurement theory. A result of fundamental importance regarding measurability of the 
electromagnetic field, proved by Bohr and Rosenfeld eighty years ago, asserts consistency of the principal limitations 
imposed by quantum theory on realizability of field measurements with formal predictions of quantum electrody- 
namics [4]. More specifically, all statistical predictions following from the formal relations between electromagnetic 
field operators can be verified, with accuracy limited only by the uncertainty principle, using an appropriately de- 
signed macroscopic measuring device. The demonstration given by Bohr and Rosenfeld refuted objections against 
the possibility of such verification, raised earlier by Landau and Peierls which were based on consideration of 
the field measurement using single test charge. It is existence of uncontrollable radiation by the test charge under 
the influence of the field being measured that led the authors of [H to the conclusion that quantum electrodynam- 
ics imposes additional limitations on the accuracy of field measurements, which turned out to be significantly more 
stringent than those following from the uncertainty principle. In particular, the use of single test charge restricts 
the accuracy of separate measurements of individual field components, in contradiction with the formal apparatus of 
quantum electrodynamics, which implies no such restriction. An important conclusion of the work [3] is that despite 
the principal impossibility to control the radiation produced by test bodies, the effect of this radiation on their motion 
can be compensated with the help of an appropriate experimental arrangement, but such compensation is possible 
only when the test bodies employed consist of sufficiently many elementary charges.^ 

Suppose now that we want to determine the electromagnetic field produced by a free electron in a given state. It 
follows from what was just said that scattering a test charge on the electron is not the best choice for this purpose. 
From the theoretical point of view, this means that using the S-matrix (say, the two-particle scattering amplitude) 
to describe the electron field generally is not adequate, as it unavoidably misses part of information about this field, 
namely that hidden by the uncontrollable radiation. Instead, since all statistical properties of the electromagnetic 
field are encoded in its operator, a complete description of the electron field can be obtained in terms of expectation 
values of the products of field operators, evaluated over the given state. Then the results of [4| guarantee that the 
predictions obtained in this way are amenable to experimental verification. One of the most important quantities of 
this sort is the expectation value of the electromagnetic field itself, called also mean, or effective, field. Thus, the 
very definition of the effective field excludes from consideration the uncontrollable radiation produced by the test 
bodies (whereas the uncontrollable radiation from the system that produces the field being measured is fully taken 
into account by the effective field formalism, Cf. Sees. Ill AllIIIAl) . 

The purpose of the present paper is to investigate infrared properties of the effective electromagnetic field in the 
case when the field-producing electron is embedded in a photon bath at finite temperature. At zero temperature, 
this problem was considered already at the dawn of quantum field theory, to determine radiative corrections to the 
field of a classical point source {e.g., atomic nucleus) The restriction to classical, that is, sufficiently heavy source 
is necessary in the conventional formulation precisely because of the presence of infrared divergences: Radiative 
corrections to the electromagnetic formfactors of charged particles vanish in the large-mass limit, which gives a formal 
reason to put the question about their divergence aside. Another reason which is often used in the literature to 
discard infrared-divergent contributions to the effective field is that at every order of the perturbation theory, such 
contributions are local, in the sense that they vanish when considered within the long-range expansion with res pec t 
to the distance from the source.'^ It is worth mentioning in this connection that despite numerous attempts |9l-[l3|. 
extension of the S-matrix formalism to finite temperatures is still an open question. Therefore, regardless of the 
fundamental reasons given above, the effective field is an indispensable means for studying temperature effects in 
quantum field theory. 

The following circumstance is crucial for the discussion of infrared singularity in the effective field. The infrared 
divergences occur because evolution of the field-producing system is considered on an infinite time interval. An 
infinite temporal extent is required already by the procedure of adiabatic switching of the interaction, which is widely 
employed in quantum field calculations, in particular, in constructing the scattering matrix. The conclusion to be 
drawn from this fact is that in the presence of massless particles, special justification is required to use the notion of 
infinitely remote past. The Bloch-Nordsieck theorem gives such justification in the case of the S-matrix, but in the 
effective field formalism this notion turns out to be physically inadequate. The point is that, as was demonstrated in 
Ref. the infrared singularity in the effective field signifies existence of a peculiar spreading of the source particle, 
which precludes preparation of a spatially localized particle state at finite times by operating with arbitrary free 
particle states in the remote past. This was shown by evaluating the effective electromagnetic field of an electron, 



^ See also Q . This result was carried over to the case of gravitational interaction by DeWitt [3 . 

^ This is actually the only reason to get rid of infrared-divergent contributions to the effective gravitational field, for the radiative 
corrections to the gravitational formfactors do not vanish in the large-mass limit. 
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regularized by means of the momentum cutoff method appropriately modified to allow factorization of the infrared 



contributions in multi-loop diagrams (called A-regularization in Ref. IJ]). The effective field was found to vanish at 
any given spatial point in the limit of removed regularization, in a way that respects the total charge conservation. 
It was also argued that the momentum cutoff can be endowed with a physical meaning as estimating the inverse 
duration of the measurement process, which made it possible to explicitly describe the electron evolution subjected 
to the irreversible spreading, and to estimate its possible observational effects. 

Except for a modification of the regularization scheme, investigation carried out in Ref. ^14"] employs the conventional 
method of calculating the effective field, based on consideration of the system on an infinite time interval. This 
investigation is therefore incomplete in view of what has been said regarding the role of the temporal extent in 
studying the infrared singularity. In addition to that, the use of an auxiliary infrared regularization raises the 
question about scheme dependence of the obtained results. At last, it is necessary to establish an exact form of 
the time-dependence of the electron state affected by the infrared singularity. Thus, we have to consider the electron 
evolution on a finite time interval, and to determine the leading large-time contributions to the electron density matrix. 
Finiteness of the time interval renders Feynman integrals infrared-convergent, removing thereby the question of scheme 
dependence, as there is no need in auxiliary infrared regularization. The infrared singularity is now contained in the 
large-time asymptotic of the effective field, or equivalently, of the effective electromagnetic current of the electron, 
and the problem is to consistently extract this asymptotic. It turns out that this requires significant modification of 
the standard calculational scheme, because restricting the consideration to a finite interval raises the issue of initial 
conditions for the electromagnetic field, which in turn enforces using the canonical Coulomb gauge instead of the 
covariant one, to avoid violation of the Gauss law in the initial state. It will be shown below how these issues are 
interrelated, and why they do not arise in the standard formulation based on the adiabatic switching of the interaction. 

The paper is organized as follows. The main tools to be used to study electron evolution in a photon bath are 
described in Sec. Ill A[ Here we also identify two time scales characterizing two essentially different stages of the electron 
evolution - the infrared thermalization corresponding to the infrared singularity, and the usual kinetic relaxation of the 
electron momentum, described by a kinetic equation. The complications introduced by finiteness of the time interval 
are discussed in detail in Sec. Ill Bl It is shown, in particular, that the standard procedure of transition from a canonical 
gauge to the Lorentz gauge cannot be accomplished in the usual way: the gauge-non-invariance of the electron density 
matrix leads to appearance of a Lorentz-non-invariant term in the Lagrangian. The resulting modification of the 
Feynman rules in the Schwinger-Keldysh method is described in Sec. IIIII The infrared thermalization is studied in 
Sec. II VI First of all, the set of diagrams contributing to the large-time asymptotic of the effective current is identified 
in Sec. IIV Al this set turns out to be different from that representing the effective current within the A-regularization 
when the initial electron state is specified at t — —oo. Factorization and summation of the infrared contributions are 
then performed in Sec. lIVBl Section |V] contains an alternative derivation of the main result, which uses specifics of 
the four-dimensional Feynman integrals to reduce the problem of extracting the large-time asymptotic to solving a 
differential equation for the electron density matrix. Some applications are given in Sec. I VII The physical meaning 
of the infrared singularity as representing an irreversible spreading of electric charges is illustrated in Sec. IVI Al by 
working out evolution of a Gaussian wave-packet. The effect of infrared thermalization on the electron diffraction in 
the classic two-slit experiment is determined in Sec. IVI Bl interpreted in terms of decoherence, and is shown to be 
detectable in a proper experimental arrangement already at room temperature. Finally, it is demonstrated in Sec. IVIII 
how interaction of the electron with non-infrared photons leads to the usual relaxation of the electron momentum. 
This is shown by using the method of Sec.|V]to obtain a differential equation for the electron momentum distribution, 
which turns out to be the usual kinetic equation. Conclusions are drawn in Sec. IVIIII The paper has an appendix 
containing discussion of the role of the non-invariant term in the Lagrangian in Lorentz gauge. 



II. GENERAL FORMULATION 



A. The model 



Consider a non-relativistic electron of mass to, interacting with virtual and real photons in equilibrium at finite 
temperature'' T <C m. We shall deal with two closely related physical quantities: the effective (mean) electromagnetic 
current of the electron, J°^, and the effective electromagnetic field produced, Fj^^ . They are defined at every space-time 



4 We use relativistic units fl — c — 1. Also, the Minkowski metric is T}f_i,i/ = diag{+l, -1,-1,-1}. 
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point X as the expectation values 

FfA^) - {FAMx)]) , Jf{x) = {J,{x)) , 

where Af^{x), J^{x) = 'ip{x)^^tp{x) are Heisenberg picture operators of the electromagnetic potential and current, 
respectively, and angular brackets denote averaging over given (mixed) state of the electron and photons. This 
definition applies to systems at finite temperatures as well as at T = 0, and according to the general rules of quantum 
theory, it yields physical quantities averaged over series of measurements in a given state. Since the electromagnetic 
field strength is linear in A^, one has F°^{x) = d^Al^{x) — d„Alf{x), where Alf{x) ~ is the effective 

electromagnetic potential. Also, the linearity of field equations implies that 

d^Ft^A = eJfix) (f) 

Before we proceed to detailed calculation of the effective quantities introduced, it is useful to express them via the elec- 
tron density matrix, which helps exposing a certain complementarity between the effective field and another important 
object - transition probability. Any state vector of the electron interacting with photons can be expanded in the direct 
products of states describing a free electron, |e), and states describing arbitrary number of free photons, \<j>), all with 
definite momenta and spin/polarizations. Consider an electron which is in a pure state jV'o) a-t the initial time instant 
to, whereas photons are in equilibrium. This system is described by the density matrix u;((/)o)|'0o)|0o)('/'o|('0o|i 
where ^(0) — e^^*/-^/ e~^*/-^, being the energy of photons in a state ji/)), and summation is over all such 
states. The probability to find the system in a state \e)\(j)) at time t > to is 

^ «;(0o ) (e I (01 C/(t , io ) I V-o ) 1 00 ) (00 1 ( V^o I C/ (io , i ) I e) 1 0) , 

where U{t,to) is the evolution operator on the interval {to,t). Being interested solely in the electron evolution, we 
sum over all final photon states. Thus, the probability of transition of the electron into a state |e) takes the form 

wiM{eMUit,toMo)\^o){M{UUito,t)\eM = (e|^.(i)|e), (2) 

4>,4>o 

where 

g{t) = ^u;(0o)(0|C/(i,io)|V'o)|0o>(0o|(V'o|C/(io,Ol0)- (3) 

Since creation of electron-positron pairs is negligible under condition T <^ m, the initial state vector |'0o)|0o) is carried 
by U {t, to) into a one-electron state. Expression ([3]) is thus nothing but the electron density matrix at time t. It is to be 
noted that as defined, this matrix takes into account all radiation effects, including those due to soft photons emitted 
by the electron (uncontrollable radiation). Indeed, summation in Eq. Q is over all possible photon states |0). Now 
consider the mean value problem, and let F be an operator built of the fermion fields only (e.g., the electromagnetic 
current). Then one has F|e)|0) = \(j))F\e), and using completeness of the products |e)|0) the effective value of F at 
time t can be written as 

F'^'it) = ^«;(0o)(0o|(V'o|(7(io,O-P^C/(i,io)|V'o)|0o) 

00 

= wi(t>o){M{MU{to,t)F\eM{(t>\{e\U{t,toMo)\4'o) 

0,0o,e 

= Y(^\g{t)F\e). (4) 

e 

It is seen that the inclusive transition probability and the effective quantities are expressed via the same object - the 
electron density matrix. An important difference is that in the first instance we deal only with diagonal elements 
of the density matrix, whereas the effective fields depend also on off-diagonal elements. From the point of view of 
the infrared problem, this means that the cancelation of infrared singularities in the S-matrix, asserted by the Bloch- 
Nordsieck theorem and its generalizations to T 7^ [l^, : implies finiteness of only diagonal elements of the 
electron density matrix in the limit to —00, t +00. In other words, consideration of the effective field gives an 
important piece of information about the system, which is not contained in the S-matrix. In fact, the existence of an 
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irreversible charge spreading, revealed by consideration of the effective electromagnetic field, is related precisely to 
the off-diagonal elements of the electron density matrix. 

In the subsequent study of the effective electromagnetic field of electron, a formally more general setting will be 
useful wherein at <o the electron is already in a mixed state described by the density matrix qq, and F is an arbitrary 
gauge invariant operator built of fermion as well as electromagnetic field operators. Equation (|4]) is then replaced by 

F'^^it) = TT,Y,wiM{M0oU{to,t)FU{t,to)\(ba) 

4>o 

= :W-iTreTr^ {Uito,t)FUit,to)goe-f'"^) 

:W=Tr4e-^^*), /3 = l/r. (5) 

where is the Hamiltonian of free photons, and Tr^, Trg, denote traces over all photon states and the single electron 
states, respectively. This expression is merely a longwinded definition F'^^(t) = {F{t)) specified to the given initial 
condition. 

Our goal is to infer physical consequences of the infrared singularity in the effective electromagnetic field of the 
electron, or equivalently, in its effective current. Accordingly, we will be interested in the large-time asymptotics of 
these quantities. There are two essentially different types of contributions which diverge for t — > oo, when treated 
within the perturbation theory. They describe different stages of the relaxation process in the system, which are 
characterized by significantly different time scales. Namely, the faster process is the infrared thermalization described 
in [l4i] within the A-regularization, which takes place in vacuum as well as at finite temperature, and which shows 
itself in a damping of the off-diagonal elements of the momentum-space density matrix of the electron. Restoring for 
a moment the ordinary units, Ref. gives the following estimate for the characteristic time of this process 



haT 



where a = e^/ {AttHc) is the fine-structure constant, and r is the characteristic length of the problem, e.g., the distance 
between the electron and the point of observation of its field, or the spacing of an interference pattern in the two-slit 
experiment (see Sec. IVIB|) . etc. Consideration of this process is the main subject of the present paper (subsequent 
calculations confirm the above estimate for ti). 

The slower process is the usual relaxation of the electron momentum, caused by its collisions with thermal photons, 
which occurs only at T 7^ 0. Its rate can be estimated by noting that the cross-section of the low-energy electron- 
photon scattering is of order (afi/mc)^, whereas the photon density is ~ {T/hc)^. Hence, the electron mean free time 
is 



T2 



This process is described by the usual kinetic equation for the electron momentum probability distribution (that is, 
the diagonal elements of the electron density matrix), and will be considered in Sec. IVIIl 
The ratio of the two time scales, 

Tl f rT\^ 

T2 



is to be considered small within the perturbation theory. Less formally, 

^^-i-(4.36rrf, 
T2 137 

where T is assumed to be expressed in kelvins, and r in centimeters, so that the ratio is small in microscopic processes 
(r ^ 10^^ cm) for all practically important temperatures. 

Below, the two stages of the electron evolution are considered separately, the reason being that, as was already 
mentioned above, they deal with different elements of the electron density matrix. 
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B. Interaction picture and gauge conditions 



To evaluate F'^^, we shall use the Schwinger-Keldysh formalism [lajlial) according to which Eq. ([5]) can be written 
in the interaction picture as 



expi I / d^yLiiy) | F{t) 



00 



(7) 



where Lj is the interaction Lagrangian, y'^-integration is along the time-contour C = Ci U C2 running from to to 
tf >t and back as shown in Fig. [1] and Tc denotes operator ordering along this contour. It is conventional to take 
the limit — > 00 in this formula, but the choice tf = t is more appropriate for our purposes. 



time 



FIG. 1: Time-integration contour in Eq. ((Tjl. 



To completely specify the scheme of calculations, one has to fix gauge invariance of the theory, which can be done 
by including a gauge-fixing term into Lagrangian density in Eq. ([T]). The amount of computational labor essentially 
depends on the gauge choice, and covariant gauges are well-known to be of great advantage over canonical ones in 
this respect, but their use requires special justification. A mere reference to gauge-independence of the effective 
electromagnetic current is not sufficient: One has to prove that the use of covariant rules gives the same results as 
the original canonical method. It turns out that Eq. (jT)) is the case where covariant techniques is not applicable, the 
failure being directly related to finiteness of the initial instant t^. To explain the point, it will be convenient to follow 
the general Faddeev-Popov method |17] of transition to covariant gauge, despite the fact that the gauge group in the 
present case is only Abelian. 

We choose to start with the Coulomb gauge (details of quantization in this gauge can be found in [Tsl. IToj) 

a,A*=0, i = l,2,3. 

This gauge admits canonical quantization, whereby the gauge condition can be considered as an operator relation. 
Expression ([7]) written in terms of the functional integral takes the form 



(8) 



expi dV(y)j F{t). 



Here d'f/'dV' is shorthand for the fermion functional integral measure, 

Jl d%lj^{x,t')<l'4)^{x,t') J]^ di;'^(x,t')dij'^{x,t'), 

and similarly for the A-field; the superscripts 1, 2 distinguish fields belonging to the corresponding branch of the time 
contour C; integration is over aU fields satisfying V'^(i) = ^j'^{t), 'ip^{t) = -0^(0, ^^(0 = ^^(i)- Next, IV^^'^), IA]'^) are 
eigenvectors of the field operators at to, e.g., Ai{to,x)\Aj) — A]{to,x)\Al), etc. Finally, the factor Be (a formally 
infinite constant) is defined by 

Be J dLjS{d,Ar) = 1, (9) 

where integration is over the gauge group, and A" is the result of action of the gauge group element lo on the field 
A^. Note that \Ai) is independent of Aq, since Aq is not a dynamical variable. Hence, integration over Aq yields the 
Coulomb law 

Ao = -A-i Jo, (10) 
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which is thus also treated as an operator relation. Now, let us define Bl (another infinite constant) by 

Bl J dLuS{d^'A%{x)-a{x)) - 1, 
with a{x) an arbitrary function of spacetime coordinates, and multiply (jS]) by unity: 



xBlJ dc.Jia^A- - a)expi d^yL{y)^ F{t). 



(11) 



A change of integration variables V''^ — > V'l '0'^ ~^ "05 A, uj ^ uj ^ gives 

F°*f(0 fdujf d^di; (V'"|^?o|V'") /" dA,(Ai'^|e-^^*|A^)Bc5(a,An 

yLiy)] Fit.x), 



xBl5(5^A^ -a)expi J d^* 



where gauge invariance of the Lagrangian L and of the operator F was taken into account. Moreover, since the 
photon state vectors are defined by transversal components of A, they are gauge-invariant, 1^^'^'^) = 1^1'^)- But this 
is not true of the fermion states IV""^'^)- Therefore, 7^ (V"^ IpoI^^), for Qo is not gauge-invariant. Therefore, 

integration over uj leads to the following expression for the effective current in the Lorentz gauge 

i?^ff(t)=K-i fdijd^ fdA, (e'^"('^)V>o|e'""('^)V'')(4'|e~'^^1A') 



xBi5(a''A^-a)expi dV(2/)j Fit), 



where a{A) = —A ^diAi. If the limit to — >■ —oo is taken, and the interaction is negligible in the remote past, the 
function (e'-^"'-^) I go le''^"^^' 0^) can be replaced by ("0"'^ |f?o|V'^)j for then the interaction can be adiabatically switched 
off in the remote past, 

e — > 0, t — > — oo, 

which kills the phase factor e"^"'^'^^. This is exactly what happens in the 5- matrix formalism, and is realized there 
through the use of in- and out-states. At finite t^, one can get rid of the phase factor by making the inverse change of 
the fermion integration variables — > e^"^"'^"^^-0, {[) — ^ e"^"'-'^^'0, but then an additional term appears in the Lagrangian 

F'^^it) = f dV-d^ / dA,{il^^\go\i^^){Al\e-f'"^\Aj) 



:BLS{d^^A^ - a)expi ^ J d^y [L{y) - e^P^i^A-^df^d^A,]^ F{t). 



Finally, independence of a allows one to average this equation with a weight exp (— i J d'^ya^{y)/2^) , ^ = const, to 
obtain 

i^-ff(t) f dV-d^ / dA, (^>o|V'')(4'|e-^''*l4'> 

(12) 



xSiCxpi dVi(2/)j F{t), 



where 

Li^L~ {d^A,f/2^ - eJ^A-i9^5,A,. (13) 
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The averaging affects only the normahzation factor denoted Jsf, as before. The Lagrangian Li generates Feynman rules 
to be used to compute the effective field in the Lorentz gauge. It is seen that in contrast to the S-matrix formalism 
(based on taking the limit Iq —oo, tf ^ +00), these rules are not Lorentz-covariant: Explicit Lorentz invariance is 
broken by the last term in Eq. (jl3p . which largely deprives Lorentz gauge of its advantage over the Coulomb gauge. 
That the non-invariant term in Li cannot be discarded is demonstrated in Appendix A where it is shown that doing 
so would violate the Gauss law already in the tree approximation.^ 
At last, it is worth mentioning that 

d^'Afix) = const. (14) 
Indeed, it follows from Eq. (IT^ that ^jf^ satisfies 

UAf{x) - (1 - ^) d.d'^Afix) = Jf{x). 

On the other hand, Eq. ([T]) is still in force, implying that d ^d'^ A'^J^ {x) = 0. 



III. PERTURBATION THEORY 



A. Feynman rules 

Perturbation expansion of Eq. ([7]) generates expressions of the form 

Tr^Tre (T^ [A^{xi)J^'{xi)A,{x2)r {X2) ■ 



which on account of comniutativity of the interaction picture operators factorize to 

Tr^ (T, [A^{xi)A,{x2) ■■■] e"^^^) Tr^ (T, [J''{xi)r {X2) ■■■]go). 

The resulting Green functions can be further expanded in the products of particle propagators using the real-time 
techniques [20,[2l[. As usual, the Tc-ordering promotes the field propagators into 2x2 matrices, e.g., 



{x~y)= iTr^ (r, [a« (x)4^) (y) 



for the electromagnetic field, and similarly for the fermion field; the matrix indices i,j take the value 1(2) for fields 
on the forward (backward) branch of the contour C. In momentum space, the electron propagator has the form 



m? — — iO ' 
D^^^\q) = 2me{-qo)6im^-q^)i^ + m), 



(15) 



where the tilde symbolizes an operation of complex conjugation with respect to which the Dirac matrices are real. 
The photon propagator reads 



1 



27rin(fc)5(fc2 



_fc2 +iO 

^I'^fc) = -27Ti[9{k„)+n(k)]6{e)d^,{k), 
D^^^\k) = -2^i[0(-feo) + n(fc)]5(fc2)d^,(ft), n{k) 



dfiu (fc) J 



W _ 1 



(16) 



^ Note that the charge conservation fpJ^ = does not entail vanishing of the non-invariant contribution: Integration by parts in Eq. I|12| l 
yields an integral of e(jQ — J^)/\~^diAi over the hyperplane j/" = to. 
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where 

rf,.(fc)=r?,.- ^"^^^'- + ^^t'^'^^'^^^^ r;^ = (l,0) (17) 

in the Coulomb gauge, whereas in the Lorentz gauge 

d,Ak) = rj,, + {^-l)^. (18) 

As usual, the interaction vertices are generated by the Lagrangian i/, and are assigned indices 1 or 2, depending on 
the branch of the contour C to which the given vertex belongs, with an extra factor (—1) for each type-2 vertex. The 
point of observation x is assigned index 1. Each propagator connects vertices of the types assigned to the propagator 
ends, and integration in the vertices is over all space and over the time interval (to, t)- Diagrammatically, the electron 
and photon propagators will be depicted by straight and wavy lines, respectively. The operators being averaged 
{Ai^{x), J^(x), etc.) will be collectively denoted by an open circle, while the interaction vertices by full circles. 

The diagrammatic representation of the effective field in the Schwinger-Keldysh technique has its specifics when 
compared to the usual S-matrix diagrammatics. As long as the single electron problem is considered, this primarily 
concerns representation of the photons present in the system. The effect of the heat bath photons is taken into account 
by the term proportional to n(fe) in Eq. As to the photons produced by the electron, their effect is also encoded 

in the internal photon lines, in contrast to the S-matrix case where photons emitted by charges (in particular, the 
uncontrollable radiation) are represented by external photon lines attached to the charged particle propagators. It 
is integration of the scattering cross-sections over momenta of these photons that cancels the infrared divergences 
due to virtual photons appearing in the loops. The difference in the graphical representation arises because in the 
S-matrix case, the photons emitted by charges are not present in the m-state, and appear only in the out-staie as a 
result of the scattering, whereas the effective field is evaluated entirely over the given m-state, and no out-state ever 
appears in the formalism. In this respect, the effective field diagrammatics bears some resemblance with that of the 
scattering cross-sections transformed using the unitarity relations. This is illustrated by Fig. HJa) which represents 
one of diagrams describing the effect of the electron-photon scattering on the electron current. The vertical broken 
line shows the cut to be made to relate this diagram with the cross-section of the Compton scattering (note that all 
photon momenta in this diagram are on the mass shell). This effect will be considered in detail in Sec. IVIIi where it 
will be shown explicitly that the expression for the electron-photon scattering cross-section is contained in diagrams 
similar to that in Fig. [2Ia) , appearing in the kinetic equation for the electron momentum probability distribution. 
Analogously, shown in Fig. [Sfb) is one of diagrams describing the effect of photon emission by the electron, that 
accompanies the scattering of heat-bath photons by the electron, etc. 




FIG. 2: Diagrams contributing to the effective electron current, (a) A fourth-order diagram describing the effect of electron- 
photon scattering, (b) A sixth-order diagram describing the effect of electron-photon scattering accompanied by a photon 
emisson by the electron. Vertical broken lines show unitarity cuts to be made to relate these diagrams to the scattering 
cross-sections. 



Because of finiteness of the time contour C, energy is not conserved in the interaction vertices. As a result, the 
usual (5-functions expressing energy conservation in the S-matrix theory become smeared, so that the vertex factor 
takes the form, in the Coulomb gauge, 

_ _ iv°t _ iv°to 

- ie^Mv) , = — (2^)3<5(3)(t;) , (19) 

m 

where v is the sum of outgoing 4-momcnta, to be referred to in what follows as the residual momentum. In the Lorentz 
gauge, the non-invariant term in Li modifies the vertex to 
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It is to be noted that the gauge invariance of F implies ^-independence of i^®*, for = 0, so that dfj,i, in the 

Lorentz gauge can be replaced by 77^1,. Moreover, regarding the analytic structure of Feynman diagrams, the Lorentz 
gauge is completely equivalent to the Coulomb gauge when F does not involve A^. In that case, there are only internal 
photon lines each of which has its ends contracted with the tensor g^{k): 

A simple calculation shows that the latter is exactly the expression (fT7|) for dajs in the Coulomb gauge. 

As a result of the smearing of the energy (5-function, the effective electromagnetic field and the effective current 
become infrared finite, so that no special regularization such as the A-regularization introduced in (l^ is needed in 
investigating the infrared effects. 

B. Pole prescriptions and multiloop Feynman integrals 




FIG. 3: A one-loop diagram contributing to the effective current. 



Regarding pole prescriptions symbolized by ±iO in Eqs. (|15l) . ((T6)) . it is worth to make the following technical 
remark. As we shall see below, evaluation of multiple Feynman integrals is often facilitated by the use of residual 
momenta as integration variables. However, this requires some care to avoid appearance of ambiguities when changing 
the order of integration. As an illustration, consider diagram in Fig. [31 Its contribution is proportional to 



d^fc f dSl f dS2 At{q2^qi-k)At{qi-q + k) ^D^^lik) 

-7 : — rzi - y^^i 



(27r)4 J (27r)4 J (27r)4 ^-(^^-iO ' m - - iO 

The imaginary infinitesimals (— iO) specify contours of integration over qf and whose independence allows inter- 
changing the order of integration. In order to go over to integration with respect to residual momenta vi, V2, we first 
change the integration variable qi — >■ Vi according to qi = q — k + Vi, and then change 52 ~^ V2 using the relation 
q2 = q + vi + V2, im which vi is treated as a complex parameter. Equation ()21|) thus takes the form 



d^k f dSi f d^V2 AtMAt{v2) ^ 



(27r)4 J (27r)4 J^^ (27r)4 m - - 1/,^ - - iO^ m-^ + ^-^^-iO^ 

where the subscript vi indicates that the position of integration contour over V2, specified by the left pole factor 
in the integrand, depends also on vi which runs the contour specified by the other pole factor. Interchanging the 
order of integration with respect to vi and V2 now is not legitimate, as it leads to ambiguity in integrating the pole 
l/[m^ — {q + vi + ^2)^] with respect to qf. But the last formula and its generalizations turn out to be useful even 
when they do not admit changing the order of integration. 
We will have more to say on this issue in Sees. IIV A[ IIVBI 



IV. INFRARED THERMALIZATION 



Upon transition to the momentum space, Eq. (|4]) with F — takes the form 

d^q d^p 



jf{x,t) = j2 



(27r)3 (27r)3' 



V2i 



/2^ 



(22) 
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where £q = \J m? + q^, po = ^g+p ^ £qi ^-nd g{t;q,q') is the electron density matrix in the momentum space 
representation; by definition, 

gaa'{t;q,q') = {e\g{t)\e'}, (23) 

where q,a are the electron momentum and spin in the state |e), and q',cr' same for |e'). This matrix is normalized 
on unity, 

Finally, the bispinor amplitudes u„{q) are also normalized on unity, UaUa — 1, and satisfy 

(jl - m)ua{q) = 0. 

The effective electromagnetic field can be found from Eq. ([1]) . Note that in the Coulomb gauge the scalar potential 
reads simply 

Using this expression and taking into account the gauge condition diAi = 0, it is not difficult to verify that the 
effective electromagnetic field satisfies Gauss law in the infinite space. 




FIG. 4: General infrared-divergent contribution to the effective current, q and p are the electron 4-momentum and 4-momentum 
transfer, respectively. 



I e ■ ' 




FIG. 5: Detailed structure of diagrams contributing to the effective current. The horizontal photon lines are supposed to be 
arbitrarily paired. 



The contributions we are interested in come from integration over small virtual and residual momenta. The results 
of imply that these contributions diverge in the limit to — > — oo, or equivalently t — oo, as In(mi) in the vacuum 
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case, and as Tt at finite temperature. Typical infrared-singular diagram is shown in Fig. |4l Its internal structure is 
further specified by Fig. [SJ where fc, k denote virtual momenta, and v, u the residual momenta; henceforth, we set 
to = 0, without restricting generality. The factor exp— i(^jt;,° -I- '^jU^)t is brought in by the current operator in 

the interaction picture, whereas factors (1 — expi?;^t), (1 — expiu^i) come from the vertex functions A^. It is easily 
verified that the product of all these factors is invariant under the following replacement 

e-i(E,«?+E,«?)t _^ 1^ At(^;) ^ At(^) ^ i '' " ^ {2^f5^^^\v) , (26) 

a fact that will be used below. 



A. Approximations 

Since we are interested in the low-energy properties of the effective current, and since we have assumed nonrelativistic 
conditions for the electron, in what follows we shall systematically neglect radiative corrections to the electron and 
photon self-energy as well as to their interaction, but only if they represent finite relative corrections to the effective 
field, that is, if they give rise to factors [1 + 0{\q\/m)] or [1 -I- 0(|p|/m)] in the leading term. Also, we shall use 
the condition T <^ m underlying our model to omit terms proportional to T/m, unless such term is leading. In this 
approximation, the momentum space density matrices at the instants t and to are related by 

Qaa' {t; q, q') = Qaa' {to\ 9, q')R{t; 9, q')e'P°* , (27) 

where the scalar factor R{t;q,q') incorporates radiative corrections (i? = 1 in the tree approximation). Indeed, the 
7-matrix structure of a A^-loop infrared-divergent diagram contributing to J*;^ is (see Fig. [5]) 



J2 q, q'W {q')l^' - - + m) • • • -1^'' - ^[^^. + i^] + m)-f^ 

y^{i-h + fi + ™)7"^ • • • - Et^' - + m)r'uAq), 

i 

where s + f = 2N. At zero temperature, the Feynman integrals diverge logarithmically, so that j^, etc. in this 
expression give rise to infrared-finite terms. For T ^ 0, the divergence is linear, and the corresponding subleading 
contribution diverges logarithmically. On dimensional grounds, its relative order is 0{T/m), so that this contribution 
can be omitted. Thus, the above expression can be replaced by 

; 9, g')^.' {(i)l^' + m) • • • 7*3/ (/ + to)7p (^ + ™)7"' • • • (^ + ™)7"= (q) 

fT.cr' 

= J2 ^?-'(*"' 9' ■ ■ ■ 29'^^2g"i • • • 2q"^u,,iq'h^u„iq), 

(T.a' 

from which Eq. (j27p follows. This relation will be derived in a more direct way in Sec. |Vl The above reasoning also 
implies a similar simplification in denominators of the electron propagators 

1 1 

etc. 



"i' - (9 - E fc. + E v^f ± iO 2<z(E h - E «0 ± iO ' 

Next, we note that graphs with a 1- vertex appearing to the left of a 2-vertex can be omitted, because they involve the 
function D^^^\q) ~ d{—qo), and therefore do not contribute at small loop momenta. In particular, all vertices on the 
incoming electron line (the right slope of diagram in Fig. [5]) must be type-1. Less trivial is the fact that all vertices 
on the outgoing electron line (the left slope of the diagram in Fig. [S]) are to be of type-2 to give rise to a nonvanishing 
contribution. To see this, imagine for a moment that the rightmost vertex on the outgoing electron line is of type-1. 
Then following the recipe formulated in Sec. IIII B| we first perform integration over uf, and find that 

du°f e"'"/* - 1 1 

— = 0, 

27ri u'j 2q{Y^ Kj + E ""j + ^ ' 

because the contour of integration can be closed in the lower half-plane of complex u° (recall that < > 0) . Therefore, 
the rightmost vertex on the left slope must be type-2, which proves the assertion, for as was shown before, vertices to 
the left of a 2-vertex must be type-2. 
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B. Factorization of infrared contributions 



To separate contributions singular in the limit i — > oo, we use the condition T ^ m to introduce a momentum 
threshold A such that T <C A <C to, which identifies the photons with < fc < A as "soft." As was already mentioned 
in Sec. IIII Al our model requires no special infrared regularization, because restriction to a finite time interval renders 
all Feynman integrals convergent at small momenta. As to ultraviolet divergences, they are supposed to be regularized 
using some conventional means, say, the Pauli-Villars technique. If, as usual, the corresponding masses are chosen 
larger than the electron mass, then the ultraviolet divergences can be isolated and subtracted without affecting 
infrared properties of the theory, and we will assume that this has been done. Using the momentum threshold, the 
perturbation series for the function q,q + p) can be written as 

oo 

R{t; q,q + p)^ Ia{p, q) (e')'^^Jv(p, q, A), (28) 

where the factor lA{p,q) is the contributions of photons with momenta fc > A, and /Ar(p, q, A) are the infrared- 
singular parts of diagrams of the type shown in Fig. |4l Introducing abridged notation Qst{k) = q^qfDl^l\k), where 
qi = q, q2 = q + p, and s, t take on values 1, 2, the functions In take the form 

r=0 ;=0 i=l ^ ' 

where N is the number of virtual photon lines, of which r (l) reside on the incoming (outgoing) electron line, while 

A 

the remaining m = N — r — I connect the two electron lines; the symbol J indicates that all loop integrals are cut 
off at fc = A; finally, the functions and incorporate electron propagators and vertex factors. They include 
sums over all permutations of vertices residing on the outgoing and incoming electron line, respectively, the factor 
r!2''/!2' accounting for redundant permutations. To put them in a form admitting factorization of multiloop diagrams, 
we proceed as explained in Sec. IIIIBI and go over to integration with respect to residual momenta. To be specific, 
consider the lowest order diagram shown in Fig. [B) Writing 



m? — — iO m? — q^ +\Q 

for the factor brought in by D^^^\ we see that the first term does not contribute (Cf. discussion at the end of 
Sec. IIV A"| . whereas the second term gives 

_ f d^fc f d^ui Atjui) f d^U2 Atiu2) , , , 

^^^-y (^^^^('V {2nr {k + u,)q2+i0j (2.)4 (.^ +.,),,+ iO + ' 

k ui 

where (ui o U2) denotes the first term with ui,U2 interchanged. We observe that the contours of integration with 
respect to uj, ""2 can be chosen so as to meet the pole prescriptions in both terms simultaneously. Namely, the itj- 
contour must go above the poles ±kq2/q2, —U2, and not intersect the M§-contour going above the poles ±fcg2/92' 
as shown in Fig. [71 Taking into account that the function Q22{k) is even, and introducing new variables wi = ui + fc, 
W2 = U2 — k yields 
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Imu 









Reu 



k{q + p) 



- < 



FIG. 7: Pole structure and contours of integration in the complex it"-plane in the integral representing F21. 



^21 



(27r)4 J (27r)4 J (27r) 



1 



At(Ml)At(u2) 

(lii + 1*2)92 + iO 
1 

(fc + ui)g2 + i0 ^ (-fc + 1*2)92 +iO 



d^fc r d^wi /■d'*W2^ At(u;i-/fc) At(«;2 + fc) 



^Q22(fc) 



(27r)4 7 (27r)4 7 (27r)4 ""^^ wig2 + iO ^292 + iO 
This consideration is readily extended to all m. The factor F^l can be written in general as {wi + ■ ■ ■ + Wk = Wk) 



y (2^ ■ ■ ■ 72^ J l\ 7^^?22(fc™+.)At(«^2.-i - WOAt 

X At(W2! + l - fci) • • • At{w2l+rn - km) ^ 



(fcm+i + W2i) 
1 



1^1(72+101^292+10 W2I+m92+iO 

where the sum is over all permutations of indices 1,2, ... ,21 + m. That this sum appears under the sign of integral 
is precisely because the imaginary infinitesimals in all its terms are of the same sign. Therefore, this sum can be 
factorized using the formula 



E 



perm 



1 



1 



1 



1 



1 



Wiq + iOW2q + iO W^q + iO 
which is easily proved by induction for all integers m. Similarly. 

Plr (^1: ■ • ■ J ^m) = 



wiq + 10 Wmq + iO ' 



(30) 



d^Wi d^W2 



■j-j- d'^kjn+i 



(2^)4 (2^)4 J 1-L (2^) 

X At(fci - W2r+l) ■ ■ ■ At{km " W2r+m) ^ 



Qll{km+i)At{km+i - W2i-l)At(-fc„i+i - W2i) 



perm 



1 



Wigi - iO VF291 - iO V1^2r-+mgi - iO 



where the sum is over all permutations of indices 1, 2, . . . , 2r + m, and factorization is accomplished using the complex 
conjugate of Eq. We thus obtain 



d^k f d^wi d'^W2 ^ At(wi-A:)At(w2 + fc) 



(27r)4 J (27r)4 (27r)' 



r022(fc) 



n 



(u'i92 + iO)(w292 + iO) 

d^W2n+i At{w2l+i - kj) 

(27r)4 W2i+zq2 + iO 
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d4fc f d^Wi A^W2 „ ,,^, At(-Wi+fc)At(-W2^fc) 

[wiqi - i0)(w2'7i - lO) 



(27r)4 y (27r)4 (27r)4 



n 



(27r)4 W2r+i9l - iO 



Finally, substitution into Eq. (|29|) gives 



/jv(p,q,A) - ^ ^ 



5 



AT 



r=0 ;=o 



9ll 922 921 

rl2^ 112' ml ~ TV! 2^ 



, 3 = .gii +522 + 2521, 



where 



Thus, 



d^k fd^wid^W2^ „, AtK-fc)At(«;2 + fc) 



(27r)4 y (27r)4 (27r)4 



(wigr + i0)(w2(7s + iO) 



q,q+p) = Ia{p, q) exp 



9 

eg 



(31) 



Using the formula 



we find 



d-^w At(w-fc) /■ dioo e-'("'«-'=«)* - 1 



(27r)4 wq + iO 



27ri Wo — fco waq^ — fcq + iO 



kq 



^VrVs 



d^k Qrsjk) 
(27r)4 {kqr){kqs) 



1 



Inserting explicit expressions (|T6|) for the photon propagator, the last expression can be put in the form 

A 

grs = VrVs 



l^^it{kqr/q°~kqjq°) 



2^ VW^rC 

itVrkqr/q" _ p-itriskqjq° 



which can be checked by inspection. Writing grs = g^-s + gJs, where g^g is the limit of g^-s for T — ^ and t fixed, the 
result of evaluation of these integrals in various limiting cases reads (assuming q <^ m, p <^ to) 



0, 



g = 



_}_P^ 

1 p^ 



In At, 



pq 



37r2 



InAi, 



mt < 1, 
mt > 1, \p\t < 1, 



(32) 



1 p2 2 p2 

— — InTi- -Tt, 

677 rrt^ „ 6T7 



1 p' 



pq 



3t7' 



\nTt- 



1 p2 



Ti < 1, 

rt> i,r|p|t/m< 1, 

■Tt, T|p|t/m>l. 



Thus, we arrive at the following expression for the large-time asymptotic of the function R{t: q,q + p) 



R{t; q,q + p)^ lip, q, t) exp - 



2ap2 



Tt 



a = — , Tt > 1. 

At7 



(33) 



(34) 
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where I(p,q,t) is the A- independent product of Ia{p,q) with the factor contributed by the term proportional to 
In Ai in Eq. p2p. Since A m, within the logarithmic accuracy this amounts to replacing In At — >■ Inmt. Conditions 
T m and Tt ^ 1 then imply that the logarithmic term is negligible in comparison with the term proportional to 
{Tt). Thus, as long as the first stage of the electron evolution is considered, one can set I{p, q,t) — 1. 

It follows from Eq. (p4l) that R{t; q, q) is infrared-finite, which in view of Eq. (|27l) means that the infrared effects do 
not change the probability distribution of the electron momentum. The vanishing of the infrared-divergent contribution 
at p = is actually a consequence of the total charge conservation, as can be seen from Eqs. ([2^ . (P7)) . It is to be 
noted also that the infrared finiteness of the diagonal elements of the electron density matrix can be viewed as a 
special case of the Bloch-Nordsieck theorem [apply this theorem to Eq. ([2])]. 



V. ALTERNATIVE DERIVATION OF EQ. (f34)l 

The exponential time dependence of the main result expressed by Eq. p4|) suggests that it can be derived from 
consideration of the density matrix evolution in infinitesimal form. We give this derivation below in order to emphasize 
the role of the approximations made, and to demonstrate a very important fact that the infrared properties of currents 
can be determined without having to sum the perturbation series explicitly. 

Let us first establish a general relation between the matrices Q{t) and go- It follows from Eq. ([3|) that in the 
momentum representation, these matrices are related by 

Q{t;q,q') = ^u;(</)o)(</)|(e|[/(t)|^o)|0o)(<^o|(^o|C/t(i)|e')l'^> 

<t>,<t>0 

= J2 «^(0o)(0|(e|C/W|eo)|</'o)(eo|eo|e[,)(</)o|(e;,|C/t(i)|e')|0), (35) 



where U{t) = U{t,0), and 



{eo\go\e'o) ^ (eo|-0o)(V'o|eo) = g{0;qo,q'„). 



For simplicity, the electron is assumed henceforth unpolarized, g^a-' ^ ^aa' , and the spin indices are omitted. By virtue 
of momentum conservation, the evolution operator has nonvanishing matrix elements only between states satisfying 

qo = q' + Q, qo = q + Q 

where Q denotes the total momentum of photons in a state and it is taken into account that this momentum is 
zero initially, Qo — 0. It follows that 

go - go = - g = p- 

Hence, the momentum space density matrices at the instants i an to = are related by 

/d^fe 
j^-:^K{t;q,p,k)Q{Q;q + k,q + p + k). (36) 

The kernel K{t\ g, p, k) can be computed using the general formulas given in Sees. [Ill IIII Al for which purpose it is 
useful to rewrite definition as 

j3„ 



g{t\ g, g') ^ j ^""^ j ^""x' ^2^:;,^^^'' - -'"''m. x')u{q')u{q)^{t, x)). (37) 

To determine the large-time behavior of g{t), one may proceed as in Sec. IIVBI prove factorization of the infrared 
contributions, and then sum the perturbation series. However, the leading term of the large-time asymptotic can be 
found more directly. As was shown in Sec. IIVB| this term is produced by interaction of the electron with equilibrium 
photons. To extract this term, we note that the property of being in equilibrium implies that the effect of such 
photons is homogeneous in time, but with one important qualification. The photon cloud surrounding the electron 
can be said to be in equilibrium at the given temperature T only when considered on time intervals S> 1/T. This is 
because at lesser time intervals (< 1/T), the quantum indeterminacy in the photon energy becomes of order of the 
photon mean energy, and in the presence of such large fiuctuations it is evidently impossible to speak about time 
homogeneity. Thus, we introduce a time tq satisfying 

Tto » 1, (38) 
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and restrict consideration to time intervals 5t > tq. This condition justifies omission of the logarithmic contributions 
due to vacuum photon-electron interaction, making thereby the photon impact on the electron evolution homogeneous 
in time. This implies that the electron density matrices at arbitrary instants t and t + 6t are related by the same 
kernel as in Eq. p6p . viz., 



/d^fc 
j^--::^K{St; q,p, k)g{t; q + k,q + p + k), 



(39) 



where K is independent of t. Whether or not this equation can be reduced to a differential equation depends on the 
structure of the kernel K{6t;q,p,k). The point is that in general, one cannot differentiate it with respect to St at 
dt = 0, because this would violate condition psp . To put it differently, if both sides of this equation are expanded 
in powers of St, the question is whether the terms 0{St^) can be neglected in comparison with the linear term. That 
this is legitimate in the case under consideration is suggested by the results of the preceding section, and is confirmed 
by the subsequent computation. Thus, we expand K{St; q,p, k) up to the first order in St, and obtain 



dg{t;q,q + p) 

dt 



d^fc 
(2^ 



where 



K{q,p,k) = —K{t;q,p,k 



K{q, p, k)g{t; q + k,q + p + k). 



d 



(40) 



(41) 



In view of the smallness of the coupling constant, ^ 1, the function K{q,p,k) can be found in the second order 
approximation. Indeed, if we use Eq. p7p to relate Q{t + St) and Q{t), where St ~ tq, then the loop expansion of 
the integral kernel K[St;q,p,k) is a power series in c^Ttq, since each loop divergences linearly, with the dominant 
contribution coming from the integration over photon momenta k ^ I/tq. Therefore, contributions of higher order in 
e are negligible, provided that tq satisfies 



e^TTo < 1. 



(42) 



It is because of the smallness of the coupling constant that this condition is consistent with Eq. (|38|) . In zeroth order 
in the coupling (free electron evolution), g{t; g, q + p) ^ e'^ *, implying that 

K{q,p,k)l^„ = ip%27r fS^^Hk). 

Diagrams representing the second-order term in K{q,p,k) are shown in Fig. [51 Using the definition (|4ip. this term 
can be written as 



dfc" 



K{q,p,k)l. = -2- I :^n{k)S{k T;,' 

Tq J 2tt [kq2)(kqi] 

+'-s^'Hk) [ ^nik')s{k-) y: ^kr^^ 

To J 2tx (k'qr)[k'qr) 



gi-rofc92/92 _ I 



~iTok'qr/q° 



The leading large-time contribution in this integral comes from k ^ I/tq, and counting powers of k in Eq. (j40p 
readily shows that this contribution is independent of tq, as expected. It is to be emphasized in this connection 
that the four-dimensionality of spacetime is essential to reach this conclusion. In a three-dimensional spacetime, for 
instance, the leading term in the function K{St; q, p, k) is quadratic in St, and transition from Eq. to Eq. PO)) 



q + p 



q+ p+k+u 



2 ^AAA/WWVl 



q + k — v\ 



(b) 



q 



(c) 



q + p 



q — v^—V2 q + p + u^+U2 



FIG. 8: Diagrams representing the second-order contribution to the integral kernel K{q,p, k) in Eq. (|40p . 
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is not legitimate. We note, finally, that under condition (|38p . fe ~ I/tq <C T can be neglected in the arguments of g 
in Eq. (|40p . so that in effect, K{q,p,k) ^ S^'^^lk). A computation similar to that performed at the end of Sec. II VB] 
then gives 



ovr m-^ 



Substituting this into Eq. pO|) . we find 

d 



m - T, — 2P 



g{t;q,q + p) , 



and therefore, 



in agreement with Eq. (p4)) . 



g, g + p) = exp ( - ^ ) e'P«*e(0; q, q + p), (43) 



VI. PHYSICAL MANIFESTATIONS OF THE INFRARED SINGULARITY 

In this section, we discuss some physical applications of the main result expressed by Eq. (^5]) . 

A. Evolution of a Gaussian wave-packet 

To illustrate the role of the irreversible charge spreading, consider an unpolarized electron prepared at io = in a 
pure state described in the momentum space by a Gaussian wavefunction 

00 (q) = [Ai^l^f^ exp{-?2q2/2} , I = const. 

This wavefunction represents the electron localized in a region of characteristic size / near the origin. The corresponding 
density matrix is go{q,q') = (l>o{q')4>o{q). Substitution of this expression in Eqs. (j43|) . ([22|) . and evaluation of the 
Gaussian integrals yields the following expression for the electron density at arbitrary instant t 

_exp|-^|, ,,.(? + _+4e,) . e.j-j (44) 

This result shows that in addition to the usual quantum spreading described by the term /m?P in /f, there is a 
spreading due to interaction of the electron with the heat-bath photons. It is easy to see that in the setting considered, 
the latter effect is dominated by the usual quantum-mechanical spreading. Indeed, for a given the minimum of 
the sum {p + t'^ /m?P) is 2t/m, which is very large compared to 49t ~ oTt/m? , since T ^ m, a ^ 1. But it is 
not difficult to give an example where the relation between the two effects is opposite. Namely, let the electron be 
prepared at to = in a pure state which is supposed to describe this electron localized at a later instant t = t va 
a region of the size I near a point Xq. Then an appropriate momentum-space amplitude in nonrelativistic quantum 
mechanics would be 



(g) = '/'o(9)e' 



where 00(9) is real; the factor e'^°'^ realizes the free electron evolution backward in time from t = r to t = 0. Taking 
(j)Q as before, the effective electron density at the instant t = t now is 



exp 



{P + 4er) 



which demonstrates that the actual size of the wave packet at t = t is Ir = {P + 40r)^/^. This simple example shows 
that the minimal uncertainty in the position of an electron evolving freely for a time t is '-^ \/Qt. Incidentally, this 



Preliminary version - October 14, 2011 



19 



fact justifies the term "irreversible spreading." It also implies that in the conventional approach based on the notion 
of infinitely remote past, the uncertainty is formally infinite, and the effective electromagnetic field is zero at any 
given spatial point. 

We note also that the minimal uncertainty is proportional to the square root of the travel time, which is characteristic 
of diffusive processes. This suggests that the irreversible spreading of the electron wave-packet can be viewed as a kind 
of Brownian motion of the electron. However, this analogy cannot be taken literally, because the irreversible spreading 
is not driven by the electron-photon collisions. This is evident from the fact that it does not lead to relaxation of the 
electron momentum - as is seen from Eq. (I43p . the momentum probability distribution g{t] q, q) is unaffected by the 
infrared singularity. Furthermore, the irreversibility of this effect shows itself also as an increase of quantum entropy 



of the electron state [IJ], which is directly related to the fact that the infrared singularity damps the off-diagonal 



elements of the electron density matrix. But in contrast to the ordinary statistics, this infrared thcrmalization takes 
place both at T ^ and T — 0, though in vacuum the process is much slower (time dependence of the effective 
electron current is exponential at finite temperature, whereas in vacuum it is a power law). These circumstances 
emphasize the specifically quantum nature of the irreversible spreading, which has no proper analog in nonrelativistic 
physics. 

How the usual relaxation of the electron momentum due to its collisions with the heat-bath photons is described 
in the present formalism will be shown in Sec. IVIIl 

Conditions considered in the last example can be experimentally realized using a magnetic lens to focus an electron 
beam, but it is more advantageous to detect the effect of infrared thcrmalization as a decoherence of the electron 
waves in the classic two-slit experiment considered in the next section. 



B. Electron diffraction 



Consider the electron diffraction in the two-slit experiment shown schematically in Fig. |9l The wavefunction of 
diffracted electrons is a superposition of two outgoing cylindric waves, which at sufficiently large distance from the 
slits has the form 

where ri (r2) is the distance between the first (second) slit and the point of observation, A is a bispinor amplitude 
independent of the spatial coordinates x,z (z being directed along the incident beam), and k = |fc| is the incident 
electron momentum. Let 2d denote the slit spacing, and L ^ d the distance between a slit and the screen at which 
we observe the interference pattern. In a vicinity of the detector {x, z L, assuming that the origin of the coordinate 





FIG. 9: Schematics of the two-slit experiment. The incident plane electron wave propagates upward; the slits are symbolized 
by small circles producing secondary electron waves. 
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system is at the screen), the wavefunction can be written as 

^\>t{x+d)^ / id _|_ ^i>t(x-d)'^ /id 



A ikL 



Then in the absence of the photon bath, the electron density is given by 

Jo(x) = i^li^o [1 + cos(>f:a;)] . 

To determine how this expression changes in the presence of thermal photons, we have to Fourier-expand the wave- 
function 



V'(x,z) = 



27rid 



+ 00 



-^oe'- / ^e 



,iq{x+d) _|_ ^iq(x-d) 



The corresponding expression for the electron density reads 

+00 

Jo{x) = ^Vd^o ^~^^^P^ [^(^'^ ~ <^)l^- (9' - <i)x\ cos{q' d) cos{qd) . 



Now, inclusion of the infrared effect of thermal photons gives for the electron density, according to Eqs. (|22|) . (|43)) . 
Jfix) 

+ 00 

= ^V'dV'o JJ exp [-er(q' - q)'^ + \d{q''^ - - i{q' - q)x\ cos(g'd) cos(gd) , 



where 8 is defined in Eq. (|44l) . and r = mL/k is the electron travel time between the slits and the detector. Integrating 
back over q,q' , we find 



J fix) = 



1 + exp 



2aTL 



) cos(>ra::) 



The exponential factor in this formula describes the decoherence caused by the infrared electron thermalization. To 
determine conditions under which this effect is appreciable, we note that yc ^ 1/r, where r is the spacing of the 
interference pattern. The interference is destroyed when the expression in the exponent becomes of order unity, or 



TL 



lO^OK/cm-eV^/^ 



where T is to be expressed in kelvins, and the electron energy e in electronvolts. The resolution of modern electron 
detectors employing magnetic lenses is a few angstrom. Therefore, for electron energy ~ 10 eV and L ~ 1 m, the effect 
is detectable already at T ^ 100 K. 



VII. RELAXATION OF THE ELECTRON MOMENTUM 



To clarify the role played by the infrared thermalization, and to better expose its distinction from the usual 
thermalization, we shall now show how interaction of the electron with non-infrared photons realizes relaxation in 
the system, i.e., how the electron momentum distribution tends to equilibrium distribution. To this end, we have to 
consider evolution of the diagonal elements of the electron density matrix, q{^\ q, q), which is described in the lowest 
order by diagrams shown in Fig. 1101 As in Sec. |Vl the leading contribution we are interested in turns out to be linear 
in time, so that dgit^ q, q)/dt can be written as a functional of g{t; q, q) at the same instant, as in Eq. (pHl) . with an 
integral kernel given exactly by diagrams in Fig. 1101 A slight change of notation in this figure is to be noted, namely, 
the electron momentum q is now off the mass shell, as it is associated with the V'-operators in Eq. p7p. symbolized in 
Fig. [TO] by open circles. On the contrary, the momentum (q — ki + ^2) is on the mass shell, as being associated with 
the external lines representing the initial density matrix. To express this fact, we write q = {(7°(fci, ^2), 9} , where 
fc2) = £q-ki+k2 +^1—^2- ^6 shall see, the leading contribution comes from integration over finite fci, fc2, such 
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that momentum q is near the mass-sheh. The two conditions {q — k\ + ^2)^ = fn^ and q^ = are clearly consistent, 
since the equation Sg^k-^+k^ — £q — + k2 has a continuum of non-trivial solutions with respect to ki ^ ^2 satisfying 
kf = k2 = 0. The latter requirement follows from the fact that all vertices on the outgoing (incoming) electron line 
are of type 2 (1), by the same reason as in Sec. IIV Al It is to be noted also that and fcj ^"^^ of the same sign, as the 
opposite would allow for a double-photon emission by a free electron (formally, the relations eq-k^+k2 = £q — k^ + k2, 
q^ = m^, kf — k2 = arc inconsistent for kik2 < 0). Hence, no singularity of the type considered presently arise at 
zero temperature, as the vacuum contribution is proportional to 9(ki)9{~k2). 

(c) „ „ (^) „ „ 

W 1, 

q — u^ — U2 

!■ 

q — k^ — u^ 

!■ 

q k I ~l~ k 2 

FIG. 10: Diagrams describing the effect of electron-photon scattering on the electron evolution. Open circles symbolize ip- 
operators in the expression p7[l for the electron density matrix. 



^AAAA/WW 1 



,1 

q + Vy +V2 



\aaaMaa/\j 1 



q — k^+Vy 





Evidently, contributions given by the disconnected diagrams in Figs.llOfc) llO( d) do not change electron momentum, 
so that Eq. (PO]) can be written as 

= / ^(9' k)s{t^ 9 + + 9) ' (45) 

where g{t,q) = g{t;q,q), C{q,k) is the connected part of the integral kernel, represented by diagrams in 
Figs. [T(ira) ITUr b). and D{q) is the disconnected part. The latter need not be calculated explicitly, since there is 
a simple relation between the two parts, following from the normalization condition (|24p . Namely, in order for this 
condition be satisfied by all solutions of Eq. (|45|) , it is necessary that 



and therefore, 

dg{t,q) f d^k 



dt J (27r)3 



[C{q, k)g{t, q + k)-C{q + fc, -k)g{t, q)] . (46) 



This is nothing but the usual kinetic equation for an electron in the photon bath, with C{q, k) playing the role of the 
probability of transition (per unit time) in which the electron goes from a state with momentum [q + k) over to a 
state with momentum q. To determine this function, we have to evaluate diagrams in Figs. fTOla) fTOFb) . It is easy to 
see that the momenta g — fci , q + k2 are off the mass shell whenever q — ki + k2 is on, and the light-like vectors ki , k2 
are nonzero. Therefore, the residual momenta can be neglected in the internal electron propagators when extracting 
the leading contribution: 



— ^ — : 2 1 stc. 



im? — {q — ki + vi) + iO ni^ — [q — ki) 

Next, integration with respect to fcj, k2 yields four terms: each diagram contributes two terms - one with fcj = +|fci|, 
^2 +1^2!, and the other with fej* = — |fci|, fcj = — 1^21- Changing fci_2 — fei,2 in the latter case, the connected 
contribution to the variation of the density matrix takes the form 



5cQ{t,q) 



d^fci r d^k2 r [du]g^ [dv]^^ (2/7r)4w(g, fci, fc2)n(fci)[l + n{k2)]g{t, q + k) 
(27r)3 y (27r)3 y [m^ - {q + vi + V2f - iO][m^ - (q - ui - U2f + iO] 

(47) 
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where 



2eq2£q+fc2|fei|2|fc2| 
xtri [t^ + J/i + m) 







x(^ + m) 



m? ~ (q ~ ki) 



r7" + 7"- 



(g-fci) 



— 2 7 +7 



1 + m 



m? — {q + 



(48) 



It is understood that q — q{ki, k2) in these formulas, the arguments of q being suppressed for brevity. To extract the 
leading large-time contribution, we rewrite Eq. (j47p as 



d^fei 

(2^ J (27r)3 7_ 



(49) 



(2/7r)£^m(g, fci + g, fc2)n(fci)[l + n(fc2)]g(t, g + fc) 

[m2 - (g(fci + ^, fc2) + + t;2)2 _ iO] [m2 - (9(fci + ^, fc2) - Ul - ^2)^ + 10] 



where ^ = (0, ^), ^ = (0, 0, ^z), and fc*^ is the root of q^ {ki, ^2) = with respect to kiz. Indeed, a shift £,z — kiz 
followed by integration over kiz removes the (5-function 5{kiz — k^z), bringing us back to Eq. (|T7)) . The leading term 
comes from integration near ^ = 0. Therefore, when extracting this term, one can set ^ = in the numerator of the 
integrand in Eq. (P^ . The function w{q, fci, fc2) given by Eq. [in which q{ki, ^2) is now on the mass shell] is then 
nothing but the probability of the scattering 

electron(g) + photon(fc2) — ^ electron(g + k) + photon(fci), 

where "electron(q)" denotes unpolarized electron with momentum q, and "photon(fc)" a photon with momentum k 
in any of the two polarization states over which summation is done for the initial as well as final photons. In the 
non-relativistic approximation. 



w{q,ki,k2) 



2m2|fci||fc2| 



1 



(fclfc2)' 



"'1"'2 



In the denominator, we expand q{ki + ^, ^2) with respect to ^ to the first order 



rn^ - {q{ki + ^, fc2) + vi + V2) 



-2q\ki,k2)[vUvUi'-^^^^) 



Introducing a new integration variable C, = ^zdq'^{ki, k2)/dkiz, and using 

g°(fci,fc2) 



Sih. ~ klJ 



dk 



Iz 



5 {q^kiM) - Sq) 



thus gives 



5cQ{t.q)=I{St) 



d^fci f d^k 



(27r)3 J (27r) 



-w{q, ki,k2)n{ki)[l + n{k2)]Q{t, q + k)5 (q"{ki, k2 



(50) 



where 



I{t) = 



2tt 



dC 



(wO + C + iO)(u" -C + iO)' 
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Integrations are done with the help of the formulas 



I'O + C + iO J 27ri J 27ri + + ( + iO 



1 _ e'C* 1 - e-<* ^ 
dC 7 7 = 27rt. 

-oo C C 

The result is I{t) = t. Substituting this into Eq. ((50)) . and dividing by St, we find the connected contribution to the 
derivative dg/dt. Comparison with Eq. (|46p now gives (we use q^{ki, ^2) = £q-ki+k^ + l^il ~ l^2|i and change notation 
fei — > p, so that k2 = p + k) 

f d^p 

C(q, ^) = y (2^^*-^' ^ + ^^1+^ + \p\-\P + k\- £q) "(P) [1 + n{p + fc)] . 
Thus, Eq. (|46)) becomes 

9) /" d'^fc /" d^p 



/d k f dp 



a< y (27r)3 J (27r)3 

X [1 + n{p + k)] g{t, q + k) - n{p + fc) [1 + n{p)] g{t, q)} . (51) 



The right-hand side of this equation is the standard form of the collision integral as representing the difference of 
"gains" and "losses" of the electron due to its collisions with photons. Substituting Eq. (|16p for the photon distribution, 
one finds that the equilibrium electron momentum distribution g{q) exp(— /Jsq). It is to be noted that no quasi- 
classic condition on the electron state has been imposed in the above derivation. In fact, the existence of irreversible 
spreading implies that the large-time evolution of the electron cannot be considered quasi-classically. 



VIII. DISCUSSION AND CONCLUSIONS 



As was explained in the introduction, the reason that makes the effective field formalism indispensable in quantum 
field theory is the principal incompleteness of the S-matrix approach regarding description of the field measurements. 
A general conclusion of our investigation is that the presence of infrared divergences in the field expectation values is 
not a sign of restricted validity of the effective field formalism, but rather an indication on inadequacy of the standard 
approach based on the assumption that the problem admits an infinite temporal extension. Specifically, we have 
shown that restricting consideration to a finite time interval makes the effective field infrared-finite, and that the 
presence of infrared divergences in the standard approach means that the radiative corrections to the classical field 
grow unboundedly with time at every order of perturbation theory. We proved also that the sum of these contributions 
is bounded, and is such that the effective electromagnetic field of an electron freely moving in a photon bath vanishes 
in the large-time limit. As the simple example given in Sec. IVI A1 demonstrates, this effect is irreversible - the electron 
coordinate variance cannot be made less than ^ aTt/m'^, if the electron travels for a time t in a photon bath at 
temperature T. This conclusion holds true also in vacuum, though the effect is much weaker in that case. Thus, 
the physical meaning of the infrared singularity in the effective electromagnetic field is the existence of irreversible 
spreading of electric charges. This conclusion was formulated in Ref. [ij] as a natural interpretation of the results 
obtained using the momentum-cutoff regularization in the standard approach. Now that essentially the same results 
have been obtained directly from the fundamental principles of quantum field theory, without having to introduce an 
auxiliary infrared regularization and to use an ad hoc definition of the effective density matrix through the effective 
field, the existence of irreversible spreading is proved to be an unequivocal consequence of quantum electrodynamics. 

It should be stressed that the irreversible spreading does not affect the scattering cross-sections themselves. As 
long as the single electron evolution in a photon bath is considered, this follows directly from Eq. ([2]) and the infrared 
finiteness of diagonal elements of the electron density matrix, proved in Sec. lIVBl But this result equally applies to any 
scattering process. The point is that the scattering amplitudes can be constructed entirely in terms of momenta (and 
polarizations) of the free particles present in the initial and final states; the standard procedure is to formally replace 
particle wave-packets by infinitely wide homogeneous beams of identical particles, erasing thereby any information 
about spatial profiles of actual particle states. In other words, it is diagonal elements of the momentum density 
matrices of particles that only matter when computing the cross-sections, and these are unaffected by the infrared 
singularity. 
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Next, some technical remarks are in order. First of all, regarding the strength of the irreversible spreading, the role 
of four-dimensionality of spacetime must be emphasized. That this factor is crucial is evident already from the relation 
of this spreading to the infrared singularities of radiative corrections, but is particularly clear from comparison with 
the usual electron-photon scattering responsible for the electron momentum relaxation considered in Sec. IVIII The 
latter is also described by the loop contributions which diverge for t — > oo (see Fig. [T0|. but by those only which are 
due to integration over finite photon momenta, and which are therefore insensitive to the spacetime dimensionality. 
Furthermore, it is the four-dimensionality of spacetime that eventually allows an infinitesimal treatment of the problem, 
given in Sec. |Vl We note in this connection that the results of Sec. |V] provide an effective tool for investigating the 
infrared problem in non-Abelian gauge theories. 

Comparison of the two stages of the electron evolution is useful also when interpreting the infrared singularity 
as a thcrmalization of the electron state. The infrared thermalization turns out to be stronger than that effected 
by the electron-photon scattering, in two respects. The infrared damping of the off-diagonal elements of the matrix 
g{t; q, q') is more complete, as it takes place for all q' ^ q, whereas the usual thermalization implies only vanishing of 
the matrix elements between states with different energy, that is, q'^ ^ q^ . The other essential difference is that the 
infrared thermalization is an 0(a)-effect, whereas the electron-photon collision effects are O(a^), which is reflected 
in the ratio of characteristic times of the two stages (see Eq. (jS])). The role of temperature in these processes is also 
quite revealing: in contrast to the usual relaxation, the infrared thermalization takes place in vacuum as well as at 
T ^ 0. All these distinctions accentuate the peculiar nature of this phenomenon which has no proper analogy in 
nonrelativistic physics. 



We thank F. A. Egorov (Russian Academy of Science), P. Pronin, A. Baurov and all participants of the theoretical 
physics seminar at the Physics Faculty of Moscow State University for their interest to our work and useful discussions. 



We saw in Sec. Ill Bl that the use of covariant techniques to compute the mean field is hindered by the fact that 
on a finite time interval, the usual procedure of adiabatic switching of the interaction is not applicable, which leads 
to appearance of a non-covariant term in the interaction Lagrangian (the last term in Eq. (jl3|) ). In this appendix, 
we demonstrate on a simple example that this complication is not a mere formality: omission of this term leads to 
violation of the Gauss law already in the tree approximation. In this approximation, the mean electromagnetic field 
is represented by the diagrams in Fig. Illl whose analytic expression is 
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Appendix A: Gauss law in Lorentz gauge 




v) 



(Al) 



where 





FIG. 11: The tree contribution to the effective electromagnetic field of electron. 
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As was explained in Sec. IIII Al we may replace here d^j^ by rj^^, so that 

1 



fc2 +10 



Substituting this into Eq. (|Aip . performing the integration over v by closing the contour of w'^-integration in the lower 
half-plane {t > 0), and using the identity p^J^{p) = we obtain 



Af{t,x)^e 



{2n) 



i\p\t 



i\p\t 



1 

2|p| \po - \p\ Po + \P 



(A2) 



The term proportional to J'^{p) is the contribution of the non-covariant term in L/. Differentiation of Eq. (jA2l) gives 



d'^Afit, x) = 0, DAf {t, x) ^ eJ^{t, x) , where J^,{t, x) 



(2^ 



Suppose now that the non-covariant term is omitted (which is equivalent to replacing 5^ — > (5^ in Eq. (jAip V Then 
the equation OA'jf — eJ^ still holds, but 



d^Af{t,x) 



d^p e 'P^sin|p|i 
^ (P) 



(27r)3 



(A3) 



It is not difficult to see that the latter equation is inconsistent with the Gauss law. For instance, consider a heavy 
particle localized in a small vicinity of the origin, so that J « 0, Jo(i, x) w 5'^'^\x). Then integration of Eq. (|A3[) and 
= eJjj, gives 



A'=**(t,a;) = 0, Af{t,x) = OU - r) 



(A4) 



where r is the distance between the charge and the observation point. Thus, at any given distance r, the Gauss law 
holds only at times t > r. In general, this law is restored only asymptotically: it is seen from Eq. (|A3|) that because 



of the factor sin \p\t^ exponentially tends to zero as t — )■ 00. 

In the tree approximation, this difficulty with the charge conservation can be overcome by imposing some special 
conditions on the photon state. However, omission of the non-invariant term in Eq. (jl3l) turns out to be much more 
harmful in higher orders of perturbation theory, namely, it leads to gauge-dependence of the effective current, which 
cannot be cured by modifying the photon state vector. 
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